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Recently, Torgerson and Mandel [Phys. Rev. Lett. 76, 3939 (1996)] have reported a disagreement 
between two schemes for measuring the phase difference of a pair of optical fields. We analyze these 
schemes and derive their associated phase-difference probability distributions, including both their 
strong and weak field limits. Our calculation confirms the main point of Torgerson and Mandel of 
the non-uniqueness of an operational definition of the phase distribution. We further discuss the 
role of postselection of data and argue that it cannot meaningfully improve the sensitivity. 
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I. INTRODUCTION 



Lack of a canonical pair for the number and phase operators, h and tp, has led to much debate in quantum mechanics 

,-H ' for many years now pj. A more recent approach to the phase question involves concentrating more on what the 

experimentalist actually does: If his goal is a precision measurement then he can perform "any" measurement followed 

Cn ■ by data analysis to extract a classical parameter p|, for instance phase shift in an interferometer. Alternatively, he 

may mentally lump together the measurement and data analysis to construct an operational phase operator for his 

specific setup. This latter approach has been championed by Mandel and his coworkers Pfl- The explicit construction 

and investigation of the properties of such phase observables can give us insight into the nature of quantum states. 

r^^. , Recently, Torgerson and Mandel Vq have compared two schemes for measuring the phase-difference between a pair 

^\ ■ of optical fields. They found that a direct scheme, where a signal is beat against a second one, and an indirect 

""^5 \ scheme, where the two signals are beat against a common local oscillator, yield different probability distributions for 

f-^ • the measured phase-difference. In particular, they found that the schemes gave radically different distributions for 

very weak signals. Torgerson and Mandel take the conflicting results as evidence of the non-uniqueness of quantum 

phase. In their analyses ambiguous data is discarded. This postselection procedure has generated some discussion in 

the literature [|]|| • We analyze these two schemes in the absence and presence of this postselection and discuss its 

interpretation. We study limits for both strong and weak signals in these schemes and give closed form expressions 

for the phase-difference probability distributions. 



II. THE NFM METHOD 

In the present section we review the operational definition of quantum phase based on the eight-port homodyne 
detector as well as its photon count statistics. 

A. NFM phase operators 

The idea of Noh, Fougeres and Mandel of operationally defined phase operators is guided by a classical analysis |f(|] 
of the eight-port homodyne interferometer |ll| shown in Fig. pi. Replacing the classical light intensities at the four 
detectors Di, (i = 3,4,5,6) by number operators hi = ajai, NFM propose [p|-|6||lC|] the phase operators 
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and 



Cm = 7W = Wtw = W = vmW 5 () 

g _ "6 ~"5 _ "65 , 2 s 

^(n 4 ~ n 3 ) 2 + (n 6 ~ h 5 ) 2 y/nl 3 + hj 5 



In the classical limit the operators Cm and Sm become JT(| the cosine and sine of the phase difference between two 
classical electromagnetic fields in the modes 1 and 2. Hence we can consider these operators to be the extension of 
this classical description of phase into the quantum domain. 

How can we calculate expectation values of a function / of these operators Cm and SW? According to the NFM- 
prescription this expectation value reads 

(f(C M ,S M ))=M Y, f( , T g= ' / T 2 W^3,n 65 ), (3) 

n43 ,wo \V n l3 + n h Vnla + nfaJ 

where W(ri4,3, n^) denotes the joint count probability for the differences 7743 = 71,4 — 71,3 and 77,65 = n 6 — n 5 at the four 
detectors. The sum in above equation is performed over all 7143 and i%65, except 77,43 — n 65 = 0. In order to avoid 
the problems associated with the definitions (Jl]) and (g) for 7743 = n^ = 0, NFM disregard such measurements and 
renormalize the joint count probability via TV = [1 — W(0, 0)] _1 . The theoretical results obtained by NFM using this 
formalism are in very nice agreement [@-f| with their experiments. 

Equation (Q) shows that the central quantity in the NFM approach is the joint count probability |12] Q 

1^(7743,77(35)= 2_^ (773,77 4 ,775,77 6 |/5out|773,774,775,«6), (4) 

where /5 ou t is the density operator of the four input modes of the interferometer after transformation through the 
beam splitters and the A/4 plate. The kct 773, 774, 775, 775) denotes the product number state corresponding to the four 
output modes and the summation J^ is performed for fixed differences 7743 and 7765. 

B. Eight-port homodyne count statistics 

We now consider the case in which the input density matrix is expressed in terms of coherent states, that is 

Pin = / d 2 /3id 2 /3 2 P(/3i,/3 2 )|/3i)i x </?i | ® |0)io io(0| ® |/3 2 > 2 2 (/3 2 | ® |0) 20 20 (0|, (5) 



where P(/3i,/3 2 ) is a two-mode Glauber-Sudarshan distribution for the fields entering ports 1 and 2. The output 
density operator just before photodetection reads 



Pout 



d 2 p 1 d 2 f3 2 P(Pi,f3 2 ) ||(/3i - &e ie )) 3 3 <|(/3i - P 2 e ie )\ 

®|i(/3i+/3 2 e l9 )) 44 (i( / 9i+/3 2 e 4e )| 
®||H0i+fte tf )) BB (§(-»& + fte w )| 

® I 2-HA - /3 2 e 4e )) 6 6 (i(-z/3i - (3 2 e w )\ , (6) 

where the phase shift 9 was introduced into the field entering port 2 for reasons which will become clear in the next 
section. 

We calculate the joint count probability for the differences 7743 and 7765 in the photocount distribution, for a given 
phase shift 9, by substituting Eq. (g|) into Eq. ([|), and obtain 

W(n i3 ,n 65 \e ie )= f d 2 ^d 2 p 2 P{f3 u P 2 )K i3 K 65 , (7) 

where the kernels -^43 and Kq§ are given by 



and 



#43 = Yl l^3||(/3i - P2e ie ))\ 2 |(n 4 |i(A +/3 2 e i9 ))| 2 
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(8) 
(9) 



These sums have been previously computed [12ll5Ulq] and we find 
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where I v {z) denotes the modified Bessel function of order v 
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C. Phase distributions from photon counts 

An eight-port homodyne detector as shown in Fig. |l| measures the photon number differences 7743 and n^ at the 
output. What we obtain are the probabilities, W(ri43, 77,65), as defined in Eq. (0), and we would like to use them to 
construct a phase distribution following the prescription given by NFM and summarized in Eqs. (|l|) and (g). This 
can be done by representing the pairs (7743,77.65) as points in a two-dimensional space where 77,43 is the x coordinate 
and 7765 is the y coordinate. Eqs. (|[) and (g) imply that a difference phase of ip corresponds to a ray in this space 
which starts at the origin and makes an angle ip with the positive x-axis. The probability assigned to ip is just the 
sum of the probabilities of the points, (7743, 77,65), which the ray passes through. 

There are two problems with this prescription. The first is that the distribution will not be smooth; it will consist 
of a number of spikes. Certain rays will intersect no (7743,7755) points yielding a value of zero for the probability of 
the corresponding angle, whereas a nearby ray will hit such a point which can lead to a nonzero probability for its 
angle. What is needed is a way to smooth the distribution and this has been provided by NFM, and we shall discuss 
it shortly. The second problem has to do with the contribution of the origin, 7743 = and 7765 = 0, to the phase 
distribution. A glance at Eqs. (p]) and (0) show that the sine and cosine are not defined at this point which implies 
that the angle is not defined either. One now has to decide what to do with the probability corresponding to this 
point when constructing the phase distribution. NFM throw out the data corresponding to this point and renormalizc 
the remaining probabilities, ^(7743, n 6 5) for 7743, n 6 5 7^ 0, accordingly. Another possibility is to associate W(0, 0) with 
a uniform distribution of the angle. The philosophy behind this approach is that because we have no information 
about how to assign W(Q, 0) to any angle , we apportion it equally among all angles. 

We follow the prescription in Sec. If A to obtain the two possible phase difference distributions 
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where 'Po(y) is the distribution in which the origin has been eliminated and Vi((p) is the distribution in which it has 
been included. In the above expressions [Eqs. (O) and (p~|)] and in all that follow one should be careful to choose 
the right branch of the function arctan according to the signs of 7743 and 77,65. 

The averaging procedure which was developed by NFM to smooth the phase distribution, and which is also used 
by Torgerson and Mandel, works in the following way. A phase shift of 9 is introduced into the field entering port 2 of 
the 8-port detector. This produces new probabilities for the photon number differences at the output which we shall 
denote by 1^(7743, 7765|e 10 ). Going back to our two-dimensional space, because the phase shift maps a phase difference 



of (p into one of ip — 8, a ray which makes an angle of ip — 6 with the x-axis is associated with a phase shift of ip. This 
gives us a difference-phase distribution for each value of the phase shift 8. 



T>o(<p\8) = (5[<p - (0 M - 0)]) = AT Yl W(n i3 ,n 65 \e ie )6 
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(15) 



and 

Vx{<p\6) = ^(0,0|e ie ) + V W(n 43 ,n 65 \e i9 )6 

"43, "65^0 

We then obtain a final phase distribution by averaging 

P(<p) = i- 1^ ' MVW) (16) 

over all of them. The 8 integrations are trivially performed, giving us 

^<P) = Y" £ ^(n 43 ,7i 65 |e J ( aretan "-/"«^), (17) 
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Vi(<p) = ^ I ^ r d8W(0,0\e t9 ) + J2 W (n 4 3,n6 5 |e l(arctan ™ 65/ " 43 - v) ) 1 , (18) 



and 



"43, "657^0 

where, as before, Vq^) does not include the point 71,43 = "65 = and Vi(ip) does. Note that in calculating Vo(<p) we 
have first averaged the unnormalized distributions and then normalized the result. The constant A/o in Eq. ( |l4| ) is 
chosen so that J Q d(pV (ip) = 1. This is the procedure adopted by NFM. However, one could just as well normalize 
the distribution for each value of 8 and then average the result. This will not give the same result as the first procedure 
because the normalization constant for each value of 9, (1 — W(0, 0|e l6, )) _1 , depends on 9. There is no obvious reason 
to prefer one procedure over the other, so we simply use the one which was chosen by NFM. 

III. SCHEMES 



Using the eight-port homodyne interferometer described in Sec. [IB , Torgerson and Mandel measure the phase 



difference between two coherent optical fields by (a) beating the two fields against each other and (b) by beating them 
against a common local oscillator. Torgerson and Mandel refer to these as the direct and indirect measurements. In this 
section we briefly describe these schemes and derive analytical expressions for the corresponding phase distributions. 

A. Direct scheme 

The direct measurement is made by beating the two input fields against each other in a single eight-port homodyne 
interferometer. The expression for the phase distribution is then simply Eq. ( |l7| ) [or Eq. (|18|)]. In what follows, we 
evaluate this expression in the limits of strong and weak fields. 

1. Strong field limit 

When one of the fields, say field 2, is strong we can use the known result for the strong local oscillator limit p3[ , 
where field 2 plays the role of the local oscillator. In this case, the phase distribution (both V and Vi) reduces to 

V(^) = ^fd 2 p 1 Jd 2 p 2 J rdrP(Jh,fo) 

x expj-[rcos<£- |/3i|cos(0 2 - <M] - [rsin^ - |/3i|sin(0 2 — (j>i)] \, (19) 

where (ij = IP^e'** , for j = 1,2. 



2. Weak fields limit 



When both fields are weak, we use the relation 
for the modified Bessel function, which allows us to write the kernels in Eq. (0) as 

( e -|(l/3i| 2 + l/32| 2 ) 
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The remaining integral in Eq. (|18|) is now trivially performed yielding 
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(21) 



i^y ^W/(0,0|e* e ) = J rf 2 /3 1 d 2 /? 2 P(/? 1 ,/3 2 )e-^ 1 l 2 +l*l 2 ). 

In this limit, the sums in Eqs. (O) and ([Lq) have only contributions from those terms in the lowest order in the 
coherent state amplitudes \/3j\. Hence only terms with 71,43 = —1,0,1 and n 65 = —1,0,1 contribute and we arrive 
finally at 



and 
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Vi(tp)~-£- I d 2 p 1 d 2 f3 2 P(l3i,fe)[l + 2\l3 1 \\p 2 \cos(<p + <j> 1 -4>2)}, 



(22) 



(23) 



up to corrections which are quartic in \@j\. Note that Eq. (g^) reproduces exactly the result obtained by Torgerson 
and Mandel M, 



mW ~27rV |fr|/|fr| + |/WI/W' 



(24) 



for input fields in coherent states |/3i) and |/3 2 ). I n this case, the phase distribution V\((p), which includes the point 
7743 = n65 = 0, reduces to 



Vi{<p)~— l+2|/3 1 ||/3 2 |cos(^ + (/.i-0 2 ) 

Z7T 



(25) 



B. Indirect scheme 



The indirect measurement scheme consists of two eight-port homodyne detectors, each performing a measurement 
of the phase distribution for one of the fields relative to a common strong local oscillator, a, as seen in Fig. g. Making 
no assumptions about the incoming state we could define the relative-phase distribution as the convolution 



r-1-K 

V{ip)= / dtpiP(,<pi,<pi-<p), (26) 

Jo 

of the joint phase distribution 

-i /*oo />OC 

p((p llVa ) = - ndn / radra^ie^.rae 4 *" I/Sine** 1 , rae 4 *") (27) 

ttJo Jo 

which is the radially integrated two-mode Q-function, in terms of the two-particle state p. We note that this reduces 
to the expected joint phase distribution 

V((p 1 ,(p 2 )=P 1 (^p 1 )P 2 (^ 2 ). (28) 

when the input fields are in the product state p — p\ ® p 2 and where the individual phase distributions Pi(<fi) have 
their obvious meaning. 

Let us now restrict our attention to the case where the input fields 1 and 2 have such a factorable state, so the 
two-mode P-function can be written P((3i,f3 2 ) = P(0i)P(P 2 ). Each phase distribution now takes the form 

Vi((Pi) = - [<PpiP(pi) f ndn 



x exp|-[r, ; cosv?i - |^|cos(0 o - 4>i)] - [nsimpi - \Pi\aia(<f> - fc)] j, (29) 

where </>o is the phase of the local oscillator and the P(j3i), i = 1, 2 are Glauber-Sudarshan distributions for the fields 
1 and 2. Evaluating the integral in Eq. (pi) we obtain the expression 

n<p) = ^J#p\p{Pi)J*p\ mK l/3i|2 - |&|2 

x ri d ri r 2 dr2e- r '- r h (2\r 1 (3i + r 2 l32e- iv \) (30) 

Jo Jo 

for the phase distribution of the indirect measurement. Again we now consider the limiting cases of this exact 
expression. 

1. Strong field limit 

In the case when the input field 2 is strong but still much weaker than the local oscillator, Eq. (|3(j) approaches 
the result Eq. (FtS) for the direct measurement very quickly. This behavior can be seen from Fig. 0, where we have 
plotted both expressions for fields 1 and 2 in coherent states and |/?2//3i| = 4. For higher values of \^ 2 \ the two curves 
coincide. 

2. Weak fields limit 

In the case when both fields are weak, the Bessel function in Eq. (Bffl) can be written as 

\z\ 2 

7 (|z|)~i + JfL 

and the integral is easily performed yielding 
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to lowest order for the phase distribution. 

We conclude this section by noting that, in contrast to the direct measurement, in this indirect measurement the 
contribution from ambiguous data is always negligible due to the presence of the strong local oscillator. 



IV. STATES WITH IDENTICAL PHASE DISTRIBUTIONS 

An unusual feature of the direct phase measurement scheme of Torgerson and Mandel (TM) is that any pair of 
two-mode input states which differ in their vacuum component, |0, 0), will have the same phase distribution pLpfl. For 
states with large photon numbers this is of little importance, but for states with very small photon numbers it has 
interesting consequences. 

Let us first briefly show that this is true. Let p(n 3 ,n4,n5,nQ\9) be the probability that D 3 detects 713 photons, D4 
detects n 4 photons, etc. (see Fig. nl), given that the phase of the input in port 2 has been shifted by 9. If the two-field 
input state is |?/>in} and U(0) describes the action of the entire apparatus, with the phase shifter in port 2 set to 6, 
then 

p{n 3 ,n 4 ,n 5j n 6 \9) = \(n 3 , n 4 , n 5 , n 6 \U(9)\ip in , 0, 0)| 2 , (32) 

where the two zeros appended to IV'in) indicate that the other two inputs are in the vacuum state. Let us express 
\ipin) as 

|Vi„)=co|0,0) + |$), (33) 

where ($|0,0) = and |co| 2 + \\&\\ 2 = 1- Because the beam splitters and the phase shifter conserve total photon 
number we have 

p(O,O,O,O|0) = |c*>| 2 , 

p(7j3,n4,n 5 ,n 6 |6>) = (n 3 , n 4 , n 5 , n 6 \U(9)\$, 0, 0)| 2 , if n 3 + n 4 + n 5 + n 6 ^ 0. (34) 



Now consider the input state 



KO = M l$> ( 35 ) 



and denote the output photon number probabilities as p'(n 3 , n 4 , 715, ne\9). (We note that this vacuum-depleted state 
\ip' in ) will be typically an entangled state when the original IV'in) is a two-mode product state.) Eq. ( p4[ ) implies that 
for n 3 + n 4 + n 5 + n 6 > 0, 

p'(n 3 ,n 4 ,n 5 ,n 6 \9) = ■jr^7^p(n 3 , n 4 , n 5 , n 6 \0). (36) 

11*11 

Because events corresponding to rij = 0, j = 1,...,4 are discarded, the probabilities p(n 3 ,n,4,n 5 , n & \9) and 
p'{n 3 , n 4 , ri5, ne\0), for n 3 + n 4 + n^ + rig > 0, determine the TM phase distributions for |^i n ) and \i/)'i n ), respec- 
tively, and they are the same up to an overall factor. When the TM phase distributions for the two states are 
computed and normalized this factor disappears. Therefore, IV'in) and \ip[ n ) have the same TM phase distribution. 

Let us see what this implies for the input state considered by Torgerson and Mandel, two weak coherent states. We 
have for |/3i|, |/?2 1 <1 

IV^in) = |/?i) <8> |/?2> - e -^ l/?l|2+l/32|2) (|0,0) + ft|l,0) + /3 2 |0, 1)). (37) 

A consequence of the result in the previous paragraph is that the state 

/ 1 X1/2 

W ta) = \ \lh\ 2 + W ) O9i|l,0) + A|0,l» (38) 

has the same TM phase distribution as |-0in) which is given by Eq. (p4|), i.e. 

-r,/ n ! d , 2cos(^ + 0i -<j> 2 ) 



2n\ |/3i|/|/3 2 | + |/3 2 |/|/3i 

This fact has been noted by NFM (jj§|. Torgerson and Mandel pointed out that the agreement between the above 
phase distribution, which agreed with their measurements, and the London distribution for l^in) is poor. Suppose we 
instead compare Eq. ( J24| ) to the London distribution for the state IV'i'n)- We find that they are identical. 

One interpretation of these results is the following. If V(<p) is taken to be the TM phase distribution of |^>; n ) one 
has the rather puzzling result that as |/?i| and |/?2| go to zero, with |/3i|/|/?2| fixed, the phase distribution does not 



change. One would expect that as \ipi n ) becomes closer to the vacuum its phase distribution would become more 
uniform. On the other hand, if V{}p) is interpreted as the phase distribution for \tp[ n ) the puzzle disappears. If \j5\\ 
and | #2 1 go to zero, with |/3i|/|/?2| fixed, the state \ip' in ) does not change and, therefore, neither should its phase 
distribution. Furthermore, its operational phase distribution (obtained via the direct measurement procedure) agrees 
with its London phase distribution. Because it discards data, the measurement procedure of Torgerson and Mandel is 
insensitive to the vacuum component of input states, which, along with the previous two observations, suggests that 
what is being measured is the phase distribution of \ip' in ) and not that of IV'in). This observation is closely related to 
one by Hradil in his discussion of the original Noh, Fougeres, and Mandel proposal ]|,[||§. 

Another objection which can be raised to the interpretation of Eq. rt24) as a phase distribution for \ipi n ) has to do 
with the connection between a phase distribution and the rotational properties of its corresponding state in phase 
space. A useful property for a phase distribution to have is providing information about how distinguishable a state is 
from a rotated version of itself. A phase distribution with this property will be able to tell us how useful a given state 
is for the measurement of phase shifts. In the case of two-mode states we would like the difference-phase distribution 
of a state \ip) to give us information about |(-0|t r D(</?o)|'0)l> which is a measure of how distinguishable \ip) is from its 
difference-phase rotated version, Uu((po)\ip). Here the operator Ud(<Po) — exp[i^(n 1 — fi 2 )] rotates mode 1 by (fio/2 
and mode 2 by —ipo/2, thereby changing the phase difference between the modes by ifo- For \?pi n ) we find 

\(fhn\U D (ipo)\iM\ * e-<l*l a+ l*l a > |l + |ft|V*°/ 2 + |/? 2 | 2 e-^°/ 2 | , (39) 

which goes to one as |/?i| and |/?2 go to zero. This shows that \ipin) and its rotated version become indistinguishable 
in this limit and \ipin) would be of little use in detecting a phase shift. This is contrary to what one sees from the TM 
phase distribution in Eq. (J24|). If \[3i\ = \p 2 \ the TM phase distributions V(<f) and V{<p+ipa) are clearly distinguishable 
if ipo is, for example, 7r/4 while according to Eq. ( |39"| ) (if \/3i\, | /3^ | <C 1) the states \ipi n ) an d fTD(vo)IV'in) are n °t- On 
the other hand, we find 



(i/L\&d(<Po)WL)\ 
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(40) 



which is independent of the size of |/3i| and |/?2| and depends only on their ratio. Eq. ( |40|) implies that if \0i\ = \{3 2 \ 
and ifo = 7r/4, then \tp' in ) and its rotated version can be distinguished which is consistent with what we expect from 
its London distribution. This again suggests interpreting the TM phase distribution in Eq. (gj) as being associated 
with |i// n ) rather than IV'in)- 

It should be noted that if all of the data is kept these questions of interpretation do not arise. In that case the 
vacuum state does contribute to the phase distribution and the phase distributions for |i/) in ) and \tp' in ) are no longer 
the same. 

V. TO POST-SELECT OR NOT? 



In Section III A we have discussed two possibilities of data analysis and have shown the phase distributions corre- 
sponding to retaining or discarding ambiguous data. In this section we briefly return to this point. 

The data analysis used in Eqs. (O) and ( |L3J ) may be compactly summarized by considering the distribution for the 
random variable Z = /(A, Y) given the distribution V{X, Y) for the random variables X and Y. This can be written 
as 



V(Z)<x f dXdY6(Z - f(X,Y))-p(X,Y). 



When the data, which in this case is a pair (A, Y), unambiguously determines Z the above formula is trivial; however, 
when the data leads to no single Z then we must make the assignment based on other considerations. We suggest 
spreading the probability evenly amongst the consistent values — uniformly across the 2ir radians. This is consistent 
with a maximum entropy assignment given no other information about the incoming state. 

In contrast the approach used in Refs. |^-|6| discards data obtained for n 43 = n 65 — 0. The excluded data does not 
appear to be very useful for informing us about phase. Thus, it is justifiable to select only the unambiguous data as has 
been done in Ref . [ 3 



However, a comparison of the behavior of the direct scheme without and with postselection for 
weak fields, Eqs. (25) and (J24|) respectively, shows that the latter has a significantly narrower distribution. We have 
argued in the previous section that one resolution of this apparent paradox is that by discarding data we are changing 
the state which we measure. Another way of looking at this issue is in terms of the sensitivity of the measurement. 



It seems that when data is discarded we achieve a greater measurement sensitivity, because the resulting phase 
distribution is narrower than would be the case if all the data were kept. How can we reconcile this with Shannon's 
information theory which teaches us that we cannot improve sensitivity by discarding information? It must be that 
the apparent difference in widths — and naively sensitivity — for these two distributions is in some sense illusory. The 
resolution to this 'paradox' is that the discarded data carries information about the overall resources used which 
would need to be factored into any meaningful measure of sensitivity. The deeper question of how we may compare 
phase distributions between different schemes, with respect to the cost of resources involved, is beyond the scope of 
this paper. However, it is clear that for states of very small photon number most measurements will yield a null result 
and will have to be thrown out. Therefore, a large number of measurements will have to be made in order to obtain 
only a few usable ones, which suggests a large expenditure of resources, and a corresponding decrease of sensitivity. 



VI. CONCLUSION 



We may generate the distribution for the relative phases of a pair of coherent states /3i and /% using the direct 
or indirect schemes described here. Our results confirm Torgerson and Mandel's claim of the non-uniqueness of 
operational phase operators; the two measurement procedures lead to different phase distributions. 

A real comparison between the schemes' sensitivity must await a careful analysis of the resources used. Notwith- 
standing this, we shall make a zeroth-order comparison here in the absence of postselection. As we noted already 
above (Section [II B 1) in the strong field limit (1 <C \02\ *C a and \/3i\ <C \02\) the two schemes, direct and indirect, go 
to the same distribution (see Fig. |3|). In the weak fields limit we have calculated, however, the situation is somewhat 
more interesting: The direct scheme reproduces the statistics of the London distribution (also called Suskind-Glogower 
or Pegg-Barnett) applied to the difference-phase 0. By contrast, the indirect scheme produces a phase distribution 
which is slightly broader by the ratio of 4/ir (~ 1.27). These facts surprised the authors: Both that one scheme is 
a realization of the London distribution in the limit where both input fields are weak coherent states, and that the 
other scheme does not perform identically. Clearly, interesting questions remain. 
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FIG. 1. Eight-port homodyne detector with the four input ports 1, 10, 2, and 20 and four output ports with detectors D3, 
D4, D5, and D§ measuring the photon count differences 7143 =714 — 713 and 7i65 =716—775. A quarter- wave plate causes a phase 
shift of 7r/2 in one arm of the interferometer. Input ports 10 and 20 remain open corresponding to a vacuum input. The input 
field 2 is phase shifted by an angle 9. A measurement of photon count differences is made for each value of 9. 
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FIG. 2. Indirect measurement scheme: A pair of eight-port homodyne detectors is used to measure the phase difference 
between the input fields 1 and 2. Each field is beat against a common strong local oscillator with coherent amplitude a. 
Each eight-port detector has four detectors measuring the photon count differences uaz and n&5. Two phase distributions are 
obtained and combined to give the phase distribution for the phase difference. 
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FIG. 3. Comparison between phase distributions for both direct (solid line) and indirect measurement (dashed line), in 
the strong field limit. In this limit, the phase distribution of the direct scheme follows from Eq. (119), whereas for the indirect 
measurement, we have used the general expression Eq. (BOJ). For both curves we have chosen the input fields 1 and 2 in coherent 
states with f3i = 1 and /?2 = —4. 
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